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ABSTRACT 



Composite failures are microscopically sequential and locally redundant. 
As a result, a composite structure reliability and its strength dependency on 
geometric size is intimately dependent on the statistics of fiber filament 
strength. A composite reliability model is needed to utilize such inherent 
materials redundancy in structural design. This investigation first establishes 
the important role of fiber diameter measurement in the characterization of 
fiber filament strength statistics and the composite reliability function, and 
second, implements the diameter measurement by laser diffraction. This 
method is automated and lends itself to industrial adoption for materials 
development, acceptance and quality control. 
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I. INTRODUCTION 



Over the past decade, the advent of composite materials has radically 
altered the field of structural design. Where the designer previously had a 
finite number of materials with which to build, an infinite number of choices 
are now available. Where design deficiencies previously could only be 
corrected by modifying the geometry of a structure, adjusting the material of 
the structure may now correct the problem. Where peculiar loading or stability 
problems previously resulted in massive, overdesigned structures, the 
technology now exists to design light and efficient ones. In short, the ability to 
design a material to meet precise specifications has unshackled the designer 
from previous conventional limitations and provided much greater flexibility. 

The capability to design the material for an application does not come 
without some costs. The ability to exploit the directional characteristics of 
composite materials also requires that many more variables be optimized 
during the design phase of a project. The manufacturing process is much more 
complex, and structural repair requires different approaches. The failure 
mechanisms of composite materials are microscopically sequential, and 
reliability, life, and strength estimates for complicated structures are difficult to 
predict. Finally, structures can be very expensive to construct, frequently 
precluding extensive destructive testing. 

Composite failures are microscopically sequential and locally redundant. 
In order to utilize microscopic redundancy of composites in structural design, a 
realistic reliability model is needed. The chain-of-bundles model [Refs. 1 and 2] 
relates fiber filament failure (which is serial) to composite failure (which is 
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locally parallel). With this model, if statistical parameters for the fiber 
filaments are accurately known, the statistical strength parameters for the 
composite can be determined. With the composite statistical strength model, 
the composite reliability dependency on the structural dimensional size and 
service stress level can be quantified. 

Traditionally, the statistical strength parameters for fibers are measured in 
terms of failure loads. This investigation first focused on the importance of 
accounting for the diameter variations in the statistical characterization of 
fibers. An analysis of the stochastic interaction of the randomness in failure 
load and in fiber diameter demonstrated the importance of fiber diameter 
measurement in the resulting composite reliability characterization. Second, 
the implementation of such diameter measurements was accomplished by an 
automated, highly accurate process involving the computer digitization and 
processing of a laser diffraction pattern from a single fiber, which is presented 
and discussed. 
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n. BACKGROUND 



Structural design of composite material structures uses the concept of 
stress, or force per area, as a method of parameters reduction. Loads can be 
measured very precisely, but accurately measuring the diameters of the fibers 
that typically make up a composite ply is difficult. Fibers frequently have 
diameters on the order of 10 |im, or approximately ten times finer than a 
human hair. 

To date, it has been common to use an average diameter for calculations. 
However, extreme value parameters such as strength, life, and ultimate stress 
are more correctly described by a statistical distribution of values. As a result, 
paradoxically, structural design and analysis are based on stress; whereas the 
materials parameters (strength) input into the analysis are based on force. This 
inconsistency may lead to erroneous and perhaps nonconservative designs. 

For example, the statistical parameters for fiber strength (mean, variability) in 
terms of stress may change when the statistical scatter in fiber diameter is 
accounted for. If the statistical distribution of the diameter of the fibers could be 
used in calculations instead of the averaged values, important parameters such 
as life, strength, and stress could be modeled much more accurately than is 
currently possible. 

A. LIFE 

Service lives of composite material structures are frequently measured in 
large time units of mean life, such as thousands or millions of years. Large 
mean lives is a consequence of large life variability, which is typically three 
decades or more. A large mean life in millions of years is needed to assure 
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high reliable life in tens of years. Since the fibers for many high performance 
composites have been synthesized only recently, no service experiences are 
available. Obviously, it is impractical to establish a statistical distribution of life 
expectancy in real time, so a suitable method of time acceleration needs to be 
established. Typically, the stress or strength of a specimen is related to the life 
of the specimen, where high strength is indicative of long life and low strength 
is indicative of short life. Phoenix and Wu [Ref. 3] showed that the level of 
stress used for life measurement is critical and that a very small change in 
stress may change the life estimate by an order of magnitude or more. 

Knowing the diameter of the fibers in a composite materials is thus crucial if 
accurate life estimates are desired. 

B. RELIABILITY 

Reliability is the probability that a structure will not fail. Usually, 
reliability is expressed in terms of failure stress or life expectancy. For small 
structures, a statistical distribution of failures can be experimentally obtained, 
an appropriate curve fit to the data, and reliability can be quantified. However, 
with large structures, statistics through destructive testing is cost prohibitive. 

One method of quantifying the reliability of large structures involves a 
sequence of probabilistic mathematical models. First, a relationship between a 
single fiber and a single composite ply (many fibers) is postulated. Similarly, a 
relationship between a single ply and a large structure can be established. In 
this way, the distribution of characteristic parameters for single fibers can be 
extrapolated to enormous structures and probabilistic reliability estimates can 
be made to useful precision. Again, because stress is dependent on the cross 
sectional area, or diameter, of the fiber, accurate diameter estimates are 
necessary. [Refs. 4-6] 
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III. DIAMETER AS A STOCHASTIC VARIABLE 



The structural design process typically uses stress, strain, and their 
combination, strength, as primary variables. Strictly speaking, strength 
parameters in terms of loads should not be used in stress analysis without 
accounting for the probabilistic interaction of the random variables. Strain is 
directly observable when tensile failure load is measured. However, because 
stress cannot be measured directly, laboratory analyses must measure force and 
area and calculate the stress using an equation (stress equals force divided by 
area). If the cross sectional area is a deterministic (rather than stochastic, or 
random) variable, the statistical strength based on stress and that based on load 
will differ only by a constant. If the cross sectional area is not deterministic, 
there will be a nonlinear relationship between stress and load distributions. 

A. OBSERVATION BY COMPUTER SIMULATION 

Appendix A is a discussion of the analytical interaction of the random 
variables failure load, cross sectional area, and failure stress. Here, computer 
simulation is used to illustrate the characteristic relationship between these 
variables. For a given set of fiber samples, both the failure loads and the cross 
sectional areas will have a statistically distributed range of values. On the 
physical grounds that a filament fails by the weakest link process, a two 
parameter Weibull distribution [Ref. 7] can be used to model the failure loads. 
The appropriate distribution function to best model the diameters of the fibers 
will vary with the manufacturing process. It may be argued that drawn fibers 
will have an upper bound diameter limited by the oriface. Other factors that 
may affect the diameter of the fibers include the rate of extrusion, the ambient 
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temperature, the viscosity of the fiber material, and the handling and loading 
of the fibers after fabrication. Rational modeling of the distribution for fiber 
diameters to account for these effects are outside the realm of this 
investigation. For illustration herein, a truncated normal distribution 
(diameter always positive) was used to represent the range of diameter values. 

A series of computer simulation experiments was conducted to simulate 
independent distributions of fiber diameter and fiber failure load using 
truncated normal and Weibull distributions respectively. Standard deviations 
for the fiber diameter distribution of 5%, 10%, and 20% were used. The stress 
was calculated for each sample and Weibull parameters were determined for 
the stress distributions [Ref. 8: pp. 426-7]. Table 1 summarizes the results of the 
simulations. The shape parameter a for stress was always less than that for 
load, and the magnitude of the difference between the two values was a 
function of the standard deviation of the fiber diameter distribution. Figure 1 
illustrates the apparent relationship observed via numerical simulation 
between fiber diameter standard deviation and the change in the shape 
parameter a. 



TABLE 1. RESULTS OF COMPUTER SIMULATION 



Std Dev. of 
Diameter, a 

0.05 



C^load 

4.979 




(onload - ^stress ®load^ 



0.0066 



0.10 



4.947 



4.855 



0.0185 



0.20 



5.033 



4.496 



0.0670 



The shape parameter of the Weibull distribution is approximately inversely 



proportional to the variability. The physical consequences of characterizing 
failure strength by load (i.e., not accounting for diameter variation) is an 
underestimate of the variability of the intrinsic strength. That is, the strength 
variability measured by load is lower than the actual strength variability 
measured by stress (accounting for diameter variation). The effect of such error 
in the shape parameter estimate on the reliability of a large structure is 
discussed in the following section. 




0.0 0.1 0.2 0.3 

Fiber Diameter Standard Deviation 



Figure 1. Effect of Fiber Diameter Standard Deviation on Shape Parameter a 

B. ESTIMATING BOUNDS FOR RELIABILITY 

Substantial amounts of research have been published on the subject of 
estimation of composite material reliability from the statistical properties of 
fiber strength. Phoenix and Smith [Ref. 5] summarized the work of Harlow and 
Phoenix [Refs. 1, 2, 9-11] and others and provided three methods of estimating 
the strength of fibrous composite materials using a chain-of-bundles model 
with local load-sharing between adjacent fibers. One of these results is used in 
this investigation; it estimated the composite statistical strength distribution 
W(x) with 
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where 



W(x) 


= min Fl k Kx), x > 0; k = 1, 2, 3, . . . 


(1) 


F[kl( x ) 


= 1 - exp { - d k ( x / P ) ka }, x > 0 


(2) 


dk 


= d k (a) = 2< k -D ( K x K 2 . . . K k . x ) « 


(3) 


Kr 


= 1 + r / 2, r = 1, 2, 3, . . . 


(4) 



[Ref. 5] 

k = 1 is the condition that when one filament fails, the entire structure fails; the ' 
representation of a single filament. Therefore, given the statistical 
representation of a single filament, i.e., F^Kx), the probability of composite 
failure can be estimated from W(x) by using Equation 1. 

The relationship between the single fiber strength distribution F^Kx) and 
the composite structure distribution W(x) is illustrated in Figure 2 (in Weibull 
probability coordinates). Three important relations are observed. 



(1) The composite distribution is no longer linear in the Weibull probability 
space (i.e., the composite failure is no longer serial). 

(2) The slope (shape parameter) of the fiber strength distribution controls the 
slope of the composite distribution at the upper tail. 

(3) Small changes in the upper tail slope leads to amplified "rotation" of the 
lower tail slope. 



The third observation is of greatest structural importance. The reliability 
of Figure 2 is normalized for the physical dimension of the interface ineffective 
length (approximately 10 Df or 50 |im). The equivalent size for a rocket motor 
case (measured in this scale) is on the order of 10 25 . That is, the portion of the 
function W(x) which controls the reliability of large structures is around 
P f = 10' 25 . There, operating at small values of ( x / (3 ) or, equivalently, with 
large structures, will necessitate accurate estimates of a if accurate distributions 
W(x) are required. 
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X 43 
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Figure 2. Sensitivity of Composite Material Reliability to Shape Parameter a 



C. EFFECT OF THE DIAMETER DISTRIBUTION ON RELIABILITY 

The stress distribution based on the results of the computer simulation for 
a 20% standard deviation in fiber diameter is plotted in Figure 3. Clearly, the 
use of a stochastic diameter variable rather than a deterministic one changes 
the value of ( x / (3 ) for a given structure size and desired reliability. If the 
standard deviation of the diameter of the fibers used to build a structure with 
fibrous composite materials was 20%, the structure would be potentially unsafe, 
because at a given structural load ( x / (3 ) the probability of failure predicted 
using the inappropriate parameter, load, is in fact lower. 

Hence, it is numerically demonstrated that ignoring the statistical 
variation in fiber diameter can lead to an erroneous shape parameter slope 
based on load, which in turn "rotates" the upper tail slope of the composite 
distribution, causing a nonconservative error in the lower tail. Because the 
lower tail is divergent, as the structure becomes larger, the error introduced by 
the ignoring the fiber diameter becomes more severe. 
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Figure 3. Reliability Eotnuate for a Fiber Diameter Standard Deviation of 20% 
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IV. DIAMETER MEASUREMENT 



Several methods of measuring the diameters of very small fibers are 
available. Each method has advantages and disadvantages. With the goal of 
transferring fiber diameter measurement from the laboratory to quantify 
assurance in materials and manufacturing, this thesis work advocates the use 
of an automated approach using laser diffraction. 

A. METHODS AVAILABLE 

Measuring the diameter of a fiber that is on the order of 10 pm requires 
careful consideration. Some methods that provide precise results may damage 
the specimen, thus trivializing the measurement process. 

1. Optical Microscope 

The optical microscope is fast, simple, inexpensive, is commonly 
available, and can produce precise results. However, it requires direct human 
manipulation and interpretation, has limited resolution, and can be fatiguing 
for the operator if large numbers of measurements are to be made. 

2. Electron Microscope 

The electron microscope produces very precise results, has excellent 
resolution, is less dependent on operator skill, and is not fatiguing for the 
operator. However, the equipment is very expensive, is not commonly 
available, and the process requires pretreating the samples with a conductive 
coating that may alter the strength of the sample. 

3. Photoconductive Cell 

Previous thesis research by Bennett [Ref. 12] demonstrated that 
locating the minimums of a diffraction pattern using a photoconductive cell 



19 



and applying these minimums to the classic Fraunhofer diffraction pattern 
theory (see Appendix B) yielded reasonable results. The method has the 
advantages of being nondestructive to the sample and of using inexpensive 
equipment. However, the precision of the measurement is dependent on the 
skill of the operator, the process is tedious, and some postprocessing calculation 
is required. 

B. AUTOMATED APPROACH USING LASER DIFFRACTION 

The diffraction theory for a slit and for a fiber are discussed in detail in 
Appendix B. These theories have been applied together with some 
inexpensive equipment to produce an automated system for measuring the 
diameter of a fiber to within 1% in less than 30 seconds. 

1. Hardware 

The following pieces of equipment were used in the automated 
diameter measuring system. 

a. MicronEye™ 

Each MicronEye™ system consists of a IS32 OpticRAM chip 
mounted behind a camera lens and connected to a computer interface device. 
Briefly, the IS32 OpticRAM chip is a conventional computer memory chip with 
its protective cover removed, exposing its rectangular array of light sensitive 
pixels. These pixels can be interfaced with a computer. The specific technical 
data for the IS32 OpticRAM, the heart of the MicronEye™ system, is contained 
in Reference 13. The horizontal resolution is better than 9 pm per pixel, and 
the chip is sensitive to light with wavelengths up to near the ultraviolet range. 
Two MicronEye™ systems were used. 

b. Macintosh™ Computer 

The Macintosh™ computer was selected for its high resolution 
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graphics capability, its ease of interface with the MicronEye™, its ease of use by 
the operator, its speed of calculation, and its ability to handle a large number of 
significant digits during calculations. 

c. Helium-Neon Laser 

A low power 1 mW Helium-Neon laser was used as a collimated 
light source. The wavelength of a Helium-Neon laser is 0.6328 pm. 

d. Digital Caliper 

A digital caliper was used to measure the distance between the 
two MicronEyes™. The caliper was able to measure distances within 0.01 mm. 

2. Software 

The Macintosh™ application CALIPER was written as part of this 
thesis research. The source code is contained in Appendix C. The 
programming language C was chosen because of its powerful graphics 
capabilities and because of its speed of execution on the Macintosh™. 

3. Equipment Setup 

Figure 4 is a schematic diagram of the system. The laser, fiber support 
stand, and MicronEyes™ were all mounted on a stiff table using rigid rails. 

The MicronEyes™ were mounted to tables on the rails that had integral 
micrometers for precise adjustment. The MicronEyes™ were connected to the 
Macintosh™ computer via interface devices. The digital caliper was mounted 
to the MicronEye™ supports. The fibers were mounted to cardboard holders 
that were in turn clamped in a support stand. The support stand was mounted 
to a table on the rails that had micrometers to allow precise adjustment in the 
directions parallel and perpendicular to the laser beam. Figure 5 shows how a 
fiber was mounted for measurement. 
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Figure 4. Schematic of an Automated Diameter Measuring System 
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Figure 5. Mounting a Fiber Sample for Measurement 
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4. Procedure 



A cardboard mounted fiber was placed in the support stand and the 
stand was adjusted to place the fiber in the center of the laser beam. The 
MicronEyes™ were then centered on diffraction pattern minimums 
symmetrically about the centerline of the laser beam. The system geometry was 
input to the CALIPER program by the operator at the beginning of program 
execution and after any movement of the components. Typically, once the 
equipment was positioned, many fibers with similar diameters could be 
measured without requiring readjustment of the equipment. 

Exposure times were adjusted until distinct minimums were located 
on the Macintosh™ screen for both MicronEyes™ (one diffraction pattern at a 
time). Figure 6 shows the desired display for one image. Again, once the 
exposure settings were determined for one fiber, they were applicable for all 
fibers with similar diameters. 




Figure 6. Centering Diffraction Minimums 

After all data had been entered and all equipment adjusted, the 
program would digitize the diffraction patterns, process them, and report the 
diameter of the fiber by three methods: using just the left-hand MicronEye™, 
using just the right-hand MicronEye™, and using an average of the two 
results. 
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5. Computer Algorithm 

If the minimums of the diffraction pattern of a fiber could be 
accurately located and measured, then the methods outlined in Appendix B 
could be applied and the diameter recovered with very high precision and 
accuracy. However, to locate these minimums exactly would require many 
exposures, and each exposure would require interaction with the operator. A 
faster and only slightly less accurate method uses only one exposure per 
MicronEye™. 

It is known that the minimum lies between points A and B on the 
diffraction pattern shown in Figure 7. The minimum is very close to half way 




center of the 

diffraction pattern 

Figure 7. Approximating the Diffraction Pattern Minimum 



between points A and B. If the diameter is calculated under the assumption 
that the minimum is located at point A, the result will be too small. If the 
diameter is calculated using point B as the minimum, the result will be too 
large. If the results are averaged, the result is very close to the actual diameter. 
Computer simulation using this method resulted in errors of less than 0.1%. 
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6. Finding the Center of the Diffraction Pattern 

The methods of Appendix B require the scattering angle 0 to be 
known. The only practical method of measuring this angle is by measuring the 
distance s from the fiber to the plane of the MicronEyes™ and the distance x 
from the centerline of the pattern to the point where the angle is desired to be 
known. Figure 8 illustrates this geometry. 



plane of the 
diffraction pattern 




Figure 8. System Geometry 

Small variations in the distance s will be negligible if s is chosen large 
with respect to x. However, small variations in x will be significant. Two 
MicronEyes™ were used to improve the resolution of the measurement of the 
distance x. Figure 9 shows how using two MicronEyes™ can locate the 
diffraction pattern center automatically, thereby reducing the error in 
measuring the distance x. 
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Storch [Ref. 14] advocated optimizing the geometry of Figure 8 to 
yield two diffraction pattern minimums per MicronEye™. This approach was 
not used for several reasons. First, an automatic diffraction pattern centering 




Xi + X2 — X + Sj + § 2 

note that x^ X / 2 + 

x 2 ^ X / 2. +5 2 

Figure 9. Locating the Center of the Diffraction Pattern 

scheme was deemed preferable to a manual one. Second, the resolution of the 
measurement of the angle 0 was improved by using a larger percentage of the 
MicronEye™ sensor area. Third, by allowing a larger distance s between the 
fiber and the MicronEye™, sensitivity to errors in measuring s could be 
reduced. A larger distance s also allowed observation of small order diffraction 
pattern minimums (n = 1, 2) further from the very intense center of the 
diffraction pattern. 

7. Summary of the Automated Approach 

The approach presented above is fast, inexpensive, automatic, simple 
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to operate, and very accurate. By properly selecting a laser based on the range of 
diameters to be measured (as discussed in Appendix B) fibers can be measured 
to within 1% or better over a range of diameters from less than 1 pm to more 
than 50pm. 
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V. RESULTS 



The results for the two sections of this thesis research are summarized 
below. 

A. DIAMETER AS A STOCHASTIC VARIABLE 

It was demonstrated that statistical parameters of composites are relatable 
to the statistical fiber strength parameter by the chain-of-bundles model [Refs. 1 
and 2]. Treating fiber diameter as a stochastic variable in reliability 
computations for composite materials consistently predicted weaker structures 
than when diameter was treated as a deterministic variable. Representing fiber 
diameter with a truncated Gaussian distribution and failure load with a two 
parameter Weibull distribution yielded a failure stress distribution that was 
accurately modelled by a two parameter Weibull distribution. Small variance 
in the shape parameter a for the distribution of a single fiber was seen to have a 
large effect on reliability predictions both for lightly loaded structures and for 
very large structures. 

B. DIAMETER MEASUREMENT 

By properly selecting a laser based on the range of diameters to be 
measured, the automatic system presented can measure fibers to within 1% or 
better over a range of diameters from less than 1 pm to more than 50 pm. 

Using an error function to directly relate Fraunhofer slit diffraction results to 
Kerker’s [Ref. 15: p. 255] method for calculating the diffraction pattern of a fiber 
produced accurate results and significantly reduced computation time and 
memory requirements. 
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VI. CONCLUSIONS 



Composite structural reliability depends strongly on the strength 
variability, or shape parameter, for the consistent single fiber filament. Very 
large composite structures with high required reliability levels are especially 
sensitive to the shape parameter of the fiber and its accurate estimation. The 
shape parameter a for failure stress depends on the stochastic interaction of 
both the failure load and the fiber diameter. Obtaining a statistical distribution 
of values for fiber diameter will necessitate an automatic, highly accurate 
measuring system. The presented system consisting of a Macintosh™ 
computer, two MicronEye™ light sensitive computer chips, and a low power 
Helium-Neon laser is an economical method of providing very accurate fiber 
diameter measurements, and hence, very accurate estimates of a. Most 
importantly, this method is automated and lends itself to industrial adoption 
for materials acceptance and quality control. 
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APPENDIX A. PROBABILISTIC INTERACTION OF MULTIPLE RANDOM 

VARIABLES 



Filament strengths are traditionally measured in terms of failure loads 
and the statistical parameters are defined in terms of load . Structural design 
and analysis are operated in terms of stress . Therefore strength parameters in 
terms of load should not be used in stress analysis without accounting for the 
probabilistic interaction of the random variables. This is acceptable only if the 
area is a deterministic variable. 

The following is an outline of how the distribution function F s (s), a 
function of the stochastic functions f D (d) and F P (p), can be obtained. 

1. VARIABLE DEFINITIONS 

The following variables are hereby identified and defined: 

P - Random (stochastic) failure load 
p - realized failure load 

A - Random area defined in terms of diameter D 

a - realized area 

D - Random diameter 

d - realized diameter 

S - Random failure stress 

s - realized failure stress 

F P - Two parameter Weibull distribution 

F d - Truncated normal distribution 

F x (x) - Cumulative Distribution Function (CDF) for random variable X 
f x (x) - Probability Density Function (PDF) for random variable X 
a, p - Parameters used to describe a Weibull statistical distribution 



30 



p, a - Parameters used to describe a Gaussian statistical distribution 



2. IMPORTANT RELATIONSHIPS 




S 


= F / A 




(A.l) 


A 


= k D 2 / 4 




(A.2) 


Fp (p) 


= Pr { P < p }, 


0 < p < OO 


(A.3) 


F D (d) 


= Pr { D < d }, 


0 < d < °° 


(A.4) 


in general 








FxM 


= Pr { X < x } 




(A.5) 


f x(x) 


s Pr { X = x } 




(A.6) 


3. GOAL 









Equation A.l and Equation A.2 relate failure stress S to failure load P and 
cross sectional area A. If P and A are deterministic variables, then calculating S 
is straightforward. If A is deterministic and P is stochastic, then 

F s (s) = F P (p) / A (A .7) 

which is again straight forward. However, if P and A are both stochastic 
variables, then the division operator has no meaning. In other words, 

F s (s) * F P (p) / F A (a) (A.8) 

because the term [ Fp(p) / F A (a) ] is undefined. Thus, the goal is to obtain an 
expression for F s (s), an unknown, in terms of known quantities. 

4. OUTLINE 

If a series of many experiments were conducted that measured both the 
failure load and the cross sectional area of a fiber, the failure stress for each fiber 
could be determined. If matching pairs of failure load P and cross sectional area 
A were plotted, the failure stress S would be represented by the slope of a line 
drawn from the origin to the point. See Figure Al. The cumulative 
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distribution function (CDF) for failure stress, F s (s), may be thought of as the 
sum of all the slopes that are less than s; this defines the domain of integration 



P 




Figure Al. Stress is the Slope of a Line Drawn From the Origin 

and is illustrated in Figure A2. Thus, F s (s) is the shaded area under the curve 
in Figure A2. 



P 




Figure Al. Illustrating the Domain of Integration 
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(A.9) 



In equation form 

f y=oo f x=ys 

F sfe) = j y=0 j x=0 fp,A< x -y) dx dy 

where f P A (p,a) is the joint probability distribution of P and A. If P and A are 
independent variables, then 

f PA (p,a) = f A (a)f P (p) (A. 10) 

Combining Equation A.9 and Equation A.10 

F s (s) = 1 J* [ f A (y) f P w] dx d y 

= j 7 f (y) [J XyS f p (x) dx] dy (A.ll) 

y=° A x=o 

However, 

f x=p f x=0 f x=p 

F p (p) = J fp (x) dx = J f p (x) dx + J fp(x) dx 

1 x=-°o x=-°° 1 x=0 

C x=p 

= J fp(x) dx (all loads presumed positive) (A.12) 

So, Equation A.ll and Equation A.12 give 

f y=oo 

F S (S) = J y=£) f A (y) F p ( ys) dy (A.13) 

Also, note that 

f A (a) = Pr { A = a } = Pr { (jc D2 / 4) = a } = Pr { D 2 = (4 a / jc) } 

= Pr { D = (4 a / k ) 0 - 5 } = f D ( (4 a / k) 0 - 5 ) (A.14) 

The fact that D > 0 eliminates negative roots. 

Choose specific distribution functions to represent D and P. 

F P (p) = 1 - exp [ - (p / (3) a ], p > 0 

(A. 15) 

where a and P are shape and scale variables, respectively (a two parameter 
Weibull distribution.) 
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= exp [ - (d - p) 2 / (2 a 2 ) ] / ( a (2 tc) 0 - 5 ), d > 0 



f D (d) 

(A. 16) 

where p is the mean and a is the standard deviation (a truncated Gaussian 
distribution). Equation A. 14 and Equation A. 16 combine to give 

f A (a) = f D ((4 a / 7t) 0 - 5 ) = exp[ - ((4a/ 7i) 0 - 5 - (i) 2 / (2 a 2 )]/ ( a (2 7t) 0 - 5 ) (A.17) 

Equation A.13, Equation A.15, and Equation A.17 can now be combined to 

give 




y=0 



1^1 -exp { - ( y s / p ) a } J dy (A.18) 

Equation A.18 can be integrated numerically for any s > 0 if the four 
constants cc, (3, p, and a are specified. The constants a and [3 come from fitting a 
two parameter Weibull distribution to the failure load data. The constants p 
and a come from fitting a truncated normal distribution to the fiber diameter 
data. 

5. SUMMARY 

The CDF for failure stress S of a fiber can be obtained from experimental 
data using stochastic representations of failure load P and fiber diameter D by 
using Equation A.18 if P and D are assumed to be independent variables. A 
truncated Gaussian distribution was used for illustration purposes to represent 
the fiber diameter, but the method can be easily extended to other statistical 
representations. 
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APPENDIX B. LASER DIFFRACTION THEORY 



This appendix discusses in detail the laser diffraction theories that were 
used during this investigation. The diffraction of a slit, the diffraction of a 
fiber, and a method of relating the two patterns are outlined and described in 
detail. Finally, the utility of these theories and their application are illustrated. 

1. FRAUNHOFER DIFFRACTION 

Characteristic diffraction patterns are produced by a sheet of laser light 
passing through a long slit when the wavelength of the laser light and the 
diameter of the slit are of the same order of magnitude. The Fraunhofer 
diffraction theory is well documented in most Physics textbooks. The following 
discussion was taken substantially from Tipler [Ref. 16: Chapter 25]. 

If the light that passes through the slit of diameter Ds is represented by N 
equally spaced point sources of equal amplitude light, then the diffraction 
pattern at any point far from the slit can be obtained by the vector summation 
of the contributions from each source. Consider Figure Bl. If the distance to 
the point P where the pattern is being calculated is sufficiently far from the slit 
so that the rays from any two sources are essentially parallel, then the path 
difference between any two sources is 

8 = d sin 0 (B.l) 

where d is the distance between the sources. The phase difference between the 
two sources is then 

Y = 2tcS/X = 2 d sin 0 / A, (B.2) 

The phase difference between any two sources that are Ds / 2 apart is 

T = 2nd sin 0 / X = 2 n [(Ds / 2) sin 0] / X = n Ds sin 0 / X (B.3) 
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Figure Bl. Modelling Fraunhofer Slit Diffraction 



If r = k , then for every source contributing to the diffraction pattern 
there is another source that is exactly 180° out of phase with it and will exactly 
cancel the contribution. The result is a minimum of the diffraction pattern, as 
illustrated in Figure B2. Letting T = k in Equation B.3 produces a more useful 
expression for a minimum of the diffraction pattern 

sin 0 = n X / Ds n = 1, 2, 3, . . . (B.4) 

where n is the number of the minimum as indexed from the center of the 
diffraction pattern. 

If the amplitude of each individual source is Ai, then the magnitude of the 
diffraction pattern at any point P far from the slit and offset from the centerline 
by an angle 0 can be calculated. Let O be the phase difference 
between the slit edges. Equation B.2 then shows that 

<D = 2 7i Ds sin 0 / A. (B.5) 
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^ 

I — Ds I 



Figure B2. Minim ums of the Fraunhofer Diffraction Pattern of a Slit 



and the total contribution to the total amplitude A from all of the individual 
sources is 
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A = A(Ds,M) = A (<t>) 




Ds/2 

A i cos ( frO) dx 
Us 

0 



= 2 Ar 



• / x \ 
sm( g ) 

O 
Ds 



X = Ds / 2 



x = 0 



= Ds Ai 



• / \ 
sm( — ) 

It 

2 



(B.6) 



Normalizing by the amplitude at the center of the diffraction pattern gives 

./Ox 1 /Ox 

sin ( — ) — cos(— ) 

Ao = Ds Ai lim — = DsAi lim — — = Ds Ai (B.7) 

<t> 0 O d>-»o 1 



A 

Ao 



sm (— ) 
2 



O 

2 



(B.8) 



Relative intensity is more useful than relative amplitude and is given by 



Io 





(B.9) 



Figure B3 is a normalized plot of the Fraunhofer diffraction pattern using 
Equation B.9. 
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Figure B3. Normalized Fraunhofer Diffraction Pattern of a Slit 



2. DIFFRACTION OF A SINGLE FIBER 



The diffraction pattern of a fiber can be calculated directly using the 
method outlined in Kerker [Ref. 15: p. 255]. If the index of refraction of the 
fiber is assumed to be perfect (m = ■=«), then the relative intensity of a fiber of 
diameter Df can be expressed by 

oo 2. 

b 0 + 2 Xb n cos ( n9 ) (B.10) 

n = 1 



J_ = 2 

Io Ko S K 



where 

Ko = 2 k / X 

9 is the scattering angle 

s is the fiber to screen distance 

X is the wavelength of the laser light 

bn = Jn(a) / Hn<2)(oc) 

a = k Df / X 

Jn(a) are Bessel functions of the first kind 
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Hn (2) (a) are Hankel functions of the second kind 
[Refs. 15 and 17] 

The bn coefficients can be simplified by some simple algebra. 
Hn®(ot) = Jn (ot) - i Yn (ot) 

where Yn (a) are Bessel functions of the second kind. 

J n < a > J>) J n (a)+iY n (a) 



H®(a) J n ( a >- i Y n(«) Jn(«) + 1 Y n (<*) 

, 2 

+ i — 



0 ») : 



J„(a)Y„(cO 



(J n (a)) 2 + (Y„(a)) 2 (J»(a)) 2 + (Y„(a))' 



(B.ll) 



— Rn + i Sn 



Rn = 



(J„fa>) 



(J n (a)) 2 + (Y„(a)) 2 



(B.12) 



(B.13) 



S n — 



J n (a) Y n (a) 



(J„(a)) 2 + (Y„(o)) 2 



Thus, the relative intensity equation can be rewritten 



Io 



— — 1 <Ro + 2 X R n cos ( n0 

oS7t L Y n = 1 / 

+ + 2 X S n cos ( n0 Jj 



(B.14) 



(B.15) 



[Ref. 14] 
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A. Minimum Number of Terms 



Figure B4 and Figure B5 illustrate the fact that as a (which is 
proportional to the fiber diameter. Equation B.10) increases, the number of 
terms n that must be computed in each summation to achieve a specified 



£ 




Figure B4. Summing Rn Within a Desired Accuracy 



£ 




Figure B5. Summing Sn Within a Desired Accuracy 
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accuracy also increases. The specified accuracy can always be achieved for the 
Rn summation with fewer terms than for the Sn summation. The number of 
terms n required to achieve an accuracy of 10' 12 for any a such that 
5 < a < 100 can be predicted from a parabolic fit of the data in Figure B5. 

Calculated = 16 + 1.3 a - 0.0013 a 2 (B.16) 

0 < (Calculated “ Cxact^ < ' — 4 (B.17) 

B. Maximum Number of Terms 

For small values of a (associated with small fibers) there exists a 
maximum number of terms that can be used in the Rn and Sn summations. 

The Bessel functions of the second kind Yn asymptotically approach -©° as a 
approaches zero. Underflow will occur at different values of a for different 
computers depending on the smallest number that can be represented. Values 
of a greater than five present no particular difficulty using double precision 
FORTRAN. 

The minimum diameter fiber that can be measured with a specified 
laser can be calculated using the wavelength of the laser. 

a = 7t Df / X (B.18) 

Df = c tX/n (B.19) 

a > 5 => ( Df ) min = (5/k)X = 1.592 X (B.20) 

Equivalently, a laser may be selected based on the minimum diameter fiber to 
be measured. 

Xmin = (ti / 5 ) Df = 0.628 Df (B.21) 

Clearly, if very small diameter fibers are to be measured and the 
magnitude of the calculations on the computer to be used appears unwieldy, 
one potential solution would be to choose a laser with a minimum 
wavelength. 
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3. COMPARISON OF THE TWO MODELS 

The Fraunhofer slit diffraction pattern and the diffraction pattern from a 
fiber are similar in shape, but are related in a nonlinear fashion. Figure B6 
shows the diffraction patterns for a slit and a fiber with the same diameters 
using Equation B.9 and Equation B.10. The distance between the minimums of 
the slit pattern and the corresponding minimums of the fiber pattern increases 
progressively as the distance from the centerline increases. 



increasing distances between minimums 




Fiber 

Slit 



Figure B6. Diffraction Patterns of a Slit and a Fiber Having the Same Diameter 



4. FIBER DIAMETER MEASUREMENT USING LASER DIFFRACTION 
Diffraction pattern minimums are easily distinguished and lend 
themselves to measurement. Measurements of the locations of the 
minimums with respect to the center of the diffraction pattern are sufficient to 
accurately compute the diameter of the fiber or slit that caused the pattern. 

One of the effects that varying Ds in Equation B.5 has is changing the 
location of the diffraction pattern minimums. Figure B 7 shows that the 
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diffraction pattern for a fiber with a diameter of 8 pm has its first minimum at 
ihe same location as the diffraction pattern for a slit with a diameter of 8.34 pm. 
This relationship between the slit and the fiber at the first diffraction pattern 
minimum can be mapped over a wide range of diameters and can be extended 
to the second minimum, etc. If the difference between the two diameters is 
expressed as a percentage of the slit diameter the relationship can be 
approximated very accurately by 

(Ds - Df) / Ds = Ki Ds *2 ‘ (B.22) 

where Ki and K2 are constants and Ds and Df are the diameters of the slit and 




Figure B7. Fiber and Slit Having First Diffraction Pattern Minimums at the 

Same Location 

the fiber. Figure B8 is plot of Equation B.22 for the first and sixth minimums of 
the diffraction patterns. Clearly, Ki and K2 are functions of the number of the 
diffraction pattern minimum as indexed from the center of the pattern. 

Equation B.22 can be solved for the fiber diameter Df 

Df = Ds [1 - Ki (n) Ds K2 <n)] (B.23) 

where n is the number of the diffraction pattern minimum. 
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Figure B8. Percentage Error in Diameter Between a Slit and a Fiber 

Note that for a given X / Ds ratio. Equation B.4 predicts a maximum 
number n of observable minimums in the diffraction pattern (because 
sin 0 < 1). Conversely, for a given node and specified X there will be a limiting 
value of Ds such that Ds > n X. A further restriction comes from Fraunhofer 
diffraction theory, which makes use of a small angle approximation 
tan 0 = sin 0 = 0 that limits 0 < 15°. For a specified X and node number n, there 
exists a minimum Ds below which the relationship between slit and fiber 
diffraction is physically invalid. In Figure B8, with n = 6 and X = 0.6328 pm, 
(Ds) min = (6)(0.6328 pm) / sin(15°) = 14.7 pm, although close correlation exists 
for Ds as small as approximately 11 pm. 

Knowing the relationship between the slit diffraction pattern and the fiber 
diffraction pattern significantly reduces the required amount of effort required 
to accurately estimate the fiber diameter. The slit diffraction pattern 
calculations are simple, and the coefficients Ki and K2 need only be calculated 
once. Table B1 summarizes the values of Ki and K2 as functions of n. 




y = 1.434e-5 * x A -0.6829 



y = 1.295e-6 * x A -0.9271 
actual data points for n = 6 
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Equation B.24, Equation B.25, and Equation B.26 summarize the polynomial 
functions Ki(n) and K2(n). 



TABLE Bl. COEFFICIENTS FOR RELATING SLIT DIFFRACTION 





TO FIBER DIFFRACTION 




n 


Ki 


K2 


1 


1.434x1 O' 5 


-0.6829 


2 


9.362 x 10- 6 


-0.7247 


3 


5.480 x 10- 6 


-0.7780 


4 


3.165 xlO- 6 


-0.8334 


5 


1.890 x 10- 6 


-0.8865 


6 


1.295 xlO- 6 


-0.9271 


(Ds - Df) / Ds 


= Ki(n)DsK2(") 


(B.24) 


Ki(n) 


= [ 17.7800 - 1.01925 n - 3.56279 n* + 


1.31954 n 3 




- 0.187208 n 4 + 9.70833e-3 n 3 ] 10^ 


(B.25) 


K2(n) 


= (-0.67780) + (2.77600e-2) n + (- 4.2691 7e-2) n* 



+ (1.130e-2) n 3 + (-1.55833e-3) n 4 + (9.0e-5) n 3 (B.26) 



where 1 < n < 6 

The most accurate results are obtained when the index n of the minimum 
being used is as small as possible. 

5. SUMMARY 

The diameter of a small fiber can be estimated to within 1 % by measuring 
the distance of the diffraction pattern minimums from the center of the pattern 
and applying Equation B.4, Equation B.23, and Equation B.24. The minimum 
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used for measurement should be as close to the center of the diffraction pattern 
as possible. The ratio X / Ds should be chosen such that 

0.001 < a/Ds) < 0.5 (B.27) 
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APPENDIX C. MACINTOSH™ COMPUTER PROGRAM 



This appendix contains the source code for the application CALIPER 
written for the Macintosh™ computer. The portions of the program that deal 
with the MicronEye™ were modified from Reference 18. Some subroutines 
were extracted essentially intact from Reference 19 and Reference 20. 
Programming ideas and techniques came from several of the references listed 
in the Bibliography. References 19 through 22 were particularly valuable 
programming references. 
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#define messagel2SRID 142 

#define messagel3SRID 143 

#define messagel4SRID 144 

#define messagel5SRID 145 

#define message!6SRID 146 
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typedef struct! 

extended value; 

Str255 text; 
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initializeO; /* initialization subroutine 

FlushEvents(everyEvent); 



InitCursorQ; /* initialize the cursor as an arrow 



<D 

o 

Tj 



05 

05 

<D 

<D 

C/5 

3 

15 

C/5 

<d 

J-i 

O 

c/5 \ 
£* 
8 c 

OJ a 
^ £ 
<D tJ 



<d 



_ w 

Ora 

2 ui s 

£ 2t qJ 

3 j2£ 

*3 ^ 

2 a,o 

C 3^3 



05 

<d 

X 

C/5 

<3 

J-l 

CD 

£ 

<d 

a—* 

05 

05 

a> 



CL) \ 

J-H * 

3^ 

T3 CD 
CD -rt 
CD 2 

o.S 

CD 

-If 

3 t: 

CD -*-* 
J-4 05 

* * 



* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* ' 

* ^ • >v 

* ©^ 

* u^r 

I 8 ’a? 

* CL, A 

* 2 o 

i it 

* £ -n 



* 

* 

* 

* 

* X3 

* *3 
> 



05 ‘v 



* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

★ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 



CD 

£ -5b 



fl3 



o 

15 



• •. • ^ 

.5 a 



a> * 
<*>1 



a> 



<n 

toO.2 <3 
«i> c 
e 'T «s 

3 

^ o ^ 



<D 

to 



J-t 

<D 

tO 

3 



to 



-'o 3-0 

£ .5 D T3 03 

£sses 



CD 



(D 



a> <d a; 

i-3 .-3 i-3 r ~. 

a— » a-* -*-* ' 4-» 

<D<D<D<D<D 
(d <d cd cd <n 



c c 3 3 3 

• • iH 

***** 



to 

£ 

15 

Q 

a» 

x 



o o « 

j2T5 2 

S.S 8 

£££ 
ctj C all i — h aD 
• >-H • •— < frl *H 



<D 

O 

*-» 

Ph 

<D 

6 

3 

05 

<D 

J-H 

'w' 

05 

bO 

?15 

SQ 

i— < U_I 



<D <D 

J-H J-. 

3 3 

05 05 

o o 

Oh Oh 
X X 
<D CD 

ri <N 

3 3 
O O 

J-H J-l 

<D <D 

S £ 
*8*2 . 



a; aj'sT 

^3’3 



05 

£ 
o 

3 3 3 . 

05 05 Jh C 

O O ^ ^ 
0,0^ ^ 
X X 

w w SfJf 

' : 2 ; 2 « 0> 
.3 .3 w w 



I 

O 

T3 

3 



ri 

2 

d 

5 

Q 

o 

3 

-*-» 

<D 

05 

"bO 

O 

3m 
Q o 

£ .2 

6 2 

<D *±3 

u >; 

II ±3 

o g 

.2 8 



•r* u/ 

QU 

<D -*-• 

X 9t 



<D 
bC 



O 

> 



57 



DisposDialog(theDialog); 

getExposureSettings(); 

reversedByteArray = (Byte *)NewPtr(256); /* 256 bytes 

fillReversedByteArray(reversedByteArray); 

SetRect(&imageContentRect,0,0 / 512,128); 
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void initExposure(theExposure) 

exposureStruct *theExposure; 
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CheckltemCscanMH^howMicronlMI/inicronl.show); 

CheckItem(scanMH,showMicron2MI / micron2.show); 
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initGeometry(theDialog); /* initialize and record the check box settings 

ShowWindow(theDialog); /* display the dialog box */ 

modalFlag = ON; /* flag to exit dialog routine */ 

while (modalFlag == ON) { /* until ok or cancel buttons are hit */ 

ModalDialog((ProcPtr)NIL,&itemHit); /* handles events in the dialog window */ 
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void initGeometry(theDialog) 

DialogPtr theDialog; 

setText(theDialog / distance_sl_ET / &distance_sl_String); 

setText(theDialog,distance_s2_ET / &distance_s2_String); 



setText(theDialog,distance_Xl_ET,&distance_Xl_String); 

setText(theDialog,distance_X2_ET / &distance_X2_String); 

setTexKtheDialog^fsEO&cfsString); 

setText(theDialog,cfxET / &cfxString); 

setText(theDialog,nodeNumberET,&:nodeNumberString); 
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void doModelessClick(itemHit) 
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void doUpdateEvent(theEvent) 
EventRecord *theEvent; 
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void doScanNow(theMode,theScanType) 

BOOLEAN theMode; /* FAST or SLOW */ 

BOOLEAN theScanType; /* SINGLE or MULTIPLE 
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displayMessageCnullSRID^ullSRID^ullSRID^ullSRID); 

InitCursor(); 

return; 



ResetAlrtStageO; 

eventMask = mDownMask + keyDownMask; 
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if (micronl.show && micron2.show) { 

if(!EventAvail(eventMask,«&theEvent)) doShortDisplay(&micronl); 
if(!EventAvail(eventMask,&theEvent)) doShortDisplay(&micron2); 
} else drawImageWindow(currentMicron); 

} while (theScanType && !EventAvail(eventMask / &theEvent) «&«& procede); 



CD 

£ 

£ 

A 

i 

£ 

O 

J-4 

u 



£ 

CD 

J-4 

1-4 

£ 

u 

£P 
• ^ > > 

•fa c 

<d 5 

^ > 
CS 

«J c- 

x S 
•gg 
Y £ 

o 

II T3 
<d £ 



£ 

O 

H£ 

£ 

• ^H 

CD 
bO 
( S 

£ 

•H 

£ 

CD 

£ 

CD 

X 



£ 

CD 

> 

CD 

0/ 

13 

T3 

Qh 

£ 

£ 

<S 

CD 

4-4 

as 

J-4 

CD 

£ 

CD 

bO 

* 



j-. 

o 

Ph 



P 

£ feO 
O 

•a 8 

£ ;rj, 

ch r 

£ CD 

u cn 



u 
CD 

£ .v 

^ CD ^ 

^ £ 

§|(3| 

i- c a> d 

*- c 60 0) 
£ « b 
5L 8 3 

fsy 

•~4 U 4^ 

'rr cd't? 
^ ^ry ^ 

fP ots o 

SPh £Ph 

r Si q> c aj 

U tn >fa cn 



4-4 

<0 



O 

CO 



* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

★ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 



O 

Oh 

CO 

4-4 

• H 

T3 

£ 

(S 

£ 

a3 

£ 

<S 

u 

CD 

£ 

O 

CD 

N 

* ^ 

13 



<s 

cd 

i 

Sd 

»£ 

4-4 

Is 

J-4 

CD 

8 

<S 

U 



o 

> 



<s 

£ 

<S 

SJ 



a> 

>% 

w 

£ 

O 



£ 

O 

• pH 

<S 

1-4 

£ 

bO-v 

cx: ^ 
£ 2 
o t 

U CD 



Pi 

w 

gb 

H W 
P cn 

So 



Q 

s 

co 

3 

£ 

v 

O 

2 

cn 

3 

£ 



£ 

(S 

X) 

o 

o 

3 

CO 

V 

c/D 

4-4 
• .-4 

pQ 

<s 

4-4 

(S 

-a 

00 

>; 



(S 

04 

o 

£ 



X> 

a, 

o 

to 



o 

o 



s s 



CD 

8 

cS 

£ 

A 

i 

(S 

J-4 

CD 

8 

(S 

U 

CD 

X 



CD 

CM 

CD 

bO 

(S 

co 

CO 

CD 



CD 

to 

(S 

co 

co 

CD 



cS 

"Op 

CO 



CO 

T3 

£ 

<S 

X> 

+ 

00 

<s 



J-4 

<S 

Ph 

o 

£ 



04 

O 

to 



£ 

O 



<s 

J-4 

£ 

bp 

v£ 

£ 

O 

u 



<D 

C 

44 

£ 

O 

J-4 

£ 

O 

•jH 

a j 

u 

£ 

o 

u 

* 



•-W 

£D 

Pi H 

£3 co 

r> 8 



4-4 QJ 

<S 

+ 



£ 
co 
Oh£ 
£ 6 



Ph aJ g 

OU" 

*£ g 

o 5 £ u pi 

UP.3 o o 



rj OO 



6 

£ 

J-4 

(S 



A 

1 

(S 

J-4 

CD 
^Ph 



CO 

£ 

£ 



Pi 

H 

CO 

S £ 

S 2 

CD 

Ph fin 
£0 

^2 o 
11 

£ 

0 

A 

1 

(S 

o3 

6 
« 

c CJ 
8 aj 
« x 



o fefeUOOS 



o 

Ph 

£ 

A 



T3 £ 



OO aj 

H H ^ *-• 
^ t: t; o 
o ^5 Ph 



£ £ 

cd o 

o A 



CD 

£ 

<s 



J-4 PO 

CD Qh 

8 n 



(S 



uu 

CD CD 

p£ p£ 



£ 

0 

A 

1 

£ 

a3 

6 
(S 

U 

HRL cu 

h Xi X 

0-) 4-4 4-4 



. . CD 

h- 4 n. 
j-4 

oO 

Ph P3 
' Ph 
+ 



CD 

J 

§<*! 



72 



theCamera->portCntrlPrms.ioCompletion = NIL; 

theCamera->portCntrlPrms.ioRefNum = theCamera->inPortIOParams.ioRefNum; 

theCamera->portCntrlPrms.CSCode =8; /* ctl code for configuration information */ 

theCamera->portCntrlPrms.csp.asncConfig = configuration; 

error = error + PBControl(&theCamera->portCntrlPrms,TRUE); /* send the ctl code & info to the port */ 
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INTEGER theSRID; 

message4 = **(StringHandle)GetString(theSRID); 

ParamText(&message4,&nullString / &nullString / &nullString); 

Alert(aIertlRID,NIL); 

ParamText(&messageO / &messagel,&message2 / &message3); 
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flipVertical(&theExposure->smooth); 
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BOOLEAN doOneExposure(theMicron / theMode) 
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exposureStruct *theExposure; 
theMicron->portCntrlPrms.CSCode = 2; 

theMicron->portCntrlPrms.ioRefN um = theMicron->inPor tIOParams.ioRefNum; 
error = PBStatus(&theMicron->portCntrlPrms,TRUE); 

if (error != noErr) OSError("bufferSizeCorrect",error,"trouble interrogating # of buffered bytes"); 
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LONGINT startTicks,tickTime; 

extended theTimeConstanbtheDelay.targetTickTime; 
theDelay = 1000.0; 

targetTickTime = theDelay * 60.0 / 1000.0; 
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for(col = 2; col < 512; col++) { 

for(row = 2; row < 129; row++) { 
bitNum = 512 * row + col; 
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negNode = col; 
negFlag = ON; 
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theMicron->outer.edge. value = fabs((edgeReference - theMicron->ou ter. pixels. value) * perPixel); 
theMicron->inner.edge.value = fabs((edgeReference - theMicron->inner.pixels.value) * perPixel); 
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ex t2Str(&formats.center / micron2.inner.center.value,&micron2.inner. center.tex 
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/**************************************************************************/ 

void calculateDiametersO ( 

displayMessageCmessage^SRID^ullSRID^ullSRID/nullSRID); 

diameter 1 = doCalculateDiameter(&micronl); 
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showl = micronl.show; 
show2 = micron2.show; 



doScanNow(SLOW,SINGLE); 
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Resource code for CALIPER © 1987. Used with RMaker. 

Written by Lieutenant Jeffrey S. Kunkel 

Naval Postgraduate School, Monterey, California 
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* 2 ;; MicronEye™ #2 button 

radioButton enabled ;; radio button dialog item, enabled 

28 6 44 129 ;; top left bottom right 



MicronEye™ #2 ;; text 

* 3 ;; exposure setting for MicronEye™ #1 

editText disabled ;; editable text dialog item, disabled 

4 139 20 187 ;; top left bottom rignt 

3.00 ;; initial setting 
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* 9 ;; label 

staticText disabled ;; static text dialog item, disabled 

122 2 138 91 ;; top left bottom right 

Inner node 2 ;; text 
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;; Inner node 2 distance to pattern center 
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18 401 50 433 ;; top left bottom right 

128 ;; CALIPER icon resource ID number 
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* 38 ;; MicronEye™ #1 & #2 diameter value 

staticText disabled ;; static text dialog item, disabled 

133 435 149 482 ;; top left bottom right 







X 




X 




X 






Q) 




jd 




,2 






X 




X 




X 






G 




G 




G 






CO 




CO 




CO 






X 




X 




X 










— -i 




V 






£ +* 




£ ~ 




£ H-» 






CD X 




0 ) X 




< 1 > 44' 






X 60 




X 60 




•- SP 






60* g 




60* G 




bo’C 






° g 




^2 g 




O g 


a; 

•H 




.2 o 
x x 




.2 o 
x ^ 




X H 






4 - o 




*4 O 




*4 O 


G 




X X 




X X 




X X 


> 




CD 4 -, 




CD 4-» 




CD 4-. 






-*-» g 




4-» ^ 




G 


^G 




o jy 


, 


cj ^a; 


, 


CJ 0 ^ 


4— » 
• 


a> 

X 


g £<x 


a; 

X 


^ x 


0 ) 

X 


X "T -h 

flj ^ X 


G 


G 


t?52£ 


3 


S 22 


3 


•“ o a> 

CO 4-» 4_» 



• N • s tv • N • s • S • N • V • V • \ • V • V • V 

• N • V • V • V • V • N • S • N • N • V • N • S • V 



60 

o 



co 

<V 

X 



G .5 
£ 

32 1 



G 

cd 



cd cd 2 
co co G 



X 

CL) 

X 

G 

G 

a; 



G 
O 

4-4 
4-1 

G 

x a> 



G box £ w 

< 1 > G 
4 _> CD 

— 1 gT 3 
bo 5:3 — 1 



O w 



O g g 

c .2 ow 

o "5 ^ a> 

s c« 

3 $ 

X _ CO 



o 3 2 B S^ 2 H 3 



C 1 h 4 > 



lu 3 

3 S 

c 3 

a; ^ x 
X <D 
MX 

<d*C ^ 

4-» CD 

gX 

g SPo E 

2^,S 

^ 2 Ow 






x> _ 

p— i G '-•h 
g O J* 

vj -J-* 

G 



a> 
60 
G 
co 

Cl,cd 



X 

CD 

X 

G 

CO 



S- 

<u X 
.ts to 
bo’C 

5 g 
2 o 

X X 

QJ .— 

■*-* VJ_ 

__ <J <1> 

CU *ZJ 'T J-. 

X « £P>< 

« to 2 2 



a> 

CO 

G 

G 

O 

V 

in 



GO 

00 

CO 

CO 

CO 



X 


X 


X 


CD 


JD 


-2 


X 


X o 


x o 


G O 


G r-H 


G r-n 


.2 ^ 


.2 m 


.52^ 


00 10 


X ON 


T3 Os 



Os 

CO 



4 -» W > 

X o 

<D r-H 
H Tf 

y oo 

X Tf 

S CO 
CO ON 



s 

=L 



© 



gs 

.a^ 

22 
CO pH 



B 

n. 



★ 



X 

H 

CJ 

* 

-i— » 

G 

4—* 

CO 



$ 

CO 

CO 



£ 

n 



ON 

<N 



v ON 
<N 



X $ 

CD CO 
~ CO 

Is 

fs 

C oo 

O (N 

H " 32 ^ 

x^o 



73 £ 
a> ® 

|s 

a> 

^ CO r 

G lo 
O CO 

cs x oo 

G pH 

* X CN C 



jD 

X 

G 

co 



X 

H 



u ^ 
X CN 
2 ^ 



H 

CO 

cd 

w 

G 

O 

G 

U 






s; 

<N 

O 



CD 

X 

co (N X 



98 



*4 ;; label 

staticText disabled ;; static text dialog item, disabled 

48 2 64 224 ;; top left bottom right 

Fiber to MicronEyes™ s (fine) ;; text 
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editText disabled ;; editable text dialog item, disabled 

72 232 88 292 ;; top left bottom right 

0.0 ;; initial value 



* 11 ;; x value 

editText disabled ;; editable text dialog item, disabled 

96 232112 292 ;; top left bottom right 

82.5 ;; initial value 
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;; alert box, one button 

/l 32 (4) ;; messageDialogRID 

;; number of items in the list 



6 



<u 



os 



73 




73 


Q) 




2 


3 




3 


03 




03 


CO 




co 




g 


'S 


cT 


<u 


cT 


£ -i-i 




C 4-1 


CD X 


• 


cd X 


a £° 
60*C 


X 

2 


.■a 6o 
60* C 


o C 




2 £ 


dial 

)ttor 

ext 


• 

2 


2 o 

73 £ 



0) X O) 
6jd£> 4 -. 60 
OS . . 03 

CO ^±5 CO 
CO CO 

a; * £<a> 
£^2 £ 



x 
a> 

VJ 4-4 C-i) -*-4 

O; X X> a; 

6oSi: bo 



a> 



os _ 

-»-» o 

co - 



OS 

co 

CXO 

■ g 



■gf 

3 £ 

.52 § 

73 JS 

£^ 

a; 4-» 

'-M 

— • <U ■ 

C'T CX 
O g-O 
u O -*-* 

•H 4L 



c 

3 -•- , o o 

.S 's .a x> 



.CO 

CO £ 

C £ 

2q c 

o O q, 
£ S~ 

tjQ o 

X a) ^ 
60q 3 
03 o 
_ co £ 
5ft £ 

2 2 d 

03 £ £ 



73 

jy 

3 

03 

co 



IS- 
Z v-G 
X--.SP 

jU 60 J-4 

u 2 £ 

r s o 
£ 

w *j O 
<U X X* 
bOOJ *- 

CO ^ 



O) 



03 



73 

3 

OS 

CO 



£ 

<u 



X 

<y 



£~ 

a; x 
.a 6c 
6 jd*E 

| 6 

.2 o 

73 



a> 
60 
03 
co 

CX$ 

2 £ 



OS 

co 4 -» O 
O) X Xi <D 



X 

CD 



60 <D 

C3 



60 



03 T" o3 03 
co ^ ^5 co 

CO CO 

CD a CX <U 

SS2 S 



73 

2 






x> 

03 

CO 




3 

03 






73 

2 






73 

2 


o 




3 


O 




•■rt 




C/5 

♦ 






3 


o 




3 


ON 




03 


ON 




w 


IN. 


73 






03 


ON 




03 


cO 




a 


CO 




v 


O 






^ s 


C 


CO 




a 


LO 


CD 


cd 


VO 




✓s 

CD 


CO 


£ 






CD 


VO) 




CD 


LO 




C 

o 


(N 

o 


>-4 


H 

u 


vo 

ON 


2 

HH 


LO 


V—/ 

CO 


C 

o 


CM 

o 


1-4 


C 

O 


o 

CN 


a 

• 


4-4 


CN 


O 


CN ^ 


o 


co £ 


CO 


cO 


4-4 
pH 4-4 


CN 


o 


CN £ 


CO 


a 


3 


CO 




03 


lo o 


o 


ON 




3 


CO 




3 


LO 


o 


JO 


vD 


<C 


4-4 

* CO 


CO? 


★ .2 


CO o 




* x> 


vO 


< 


* o 


CO 


U 



103 



* 3 ;; message static text item 

staticText disabled ;; static text dialog item, disabled 

3 50 96 307 ;; top left bottom right 

A 0 A 1 ;; message text 



iconltem disabled ;; icon item, disabled 

3 9 35 41 ;; top left bottom right 

0 ;; stop icon 
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iconltem disabled ;; icon item, disabled 

3 9 35 41 ;; top left bottom right 

0 ;; stop icon 
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;; error message dialog 

,9999 ;; resource ID 

A Message... ;; message (not used) 

87 67 285 478 ;; top left bottom right 

Invisible NoGoAway ;; box status 
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Type STR ;; String Resources 

,128 ;; nullSRID 

\20 ;; ASCII code for blank space 
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One of the exposure settings contains an ++ 
illegal character. The exposure settings must ++ 
all be between zero and 1000. 
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At least one of the inputs ++ 
is still either negative or ++ 
too large, or it contains an ++ 
illegal character. 



,143 ;; messagel3SRID 

At least one of the inputs ++ 
is still wrong. 
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